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We consider the Burgers-type system studied by Foursov, 

wt = Wxx + ^wwx + (2 - 'ia)zzx, 



l^' zt = {l-2a)zxx- '^azwx + {'i-'&a)wzx- {'i + Sa)w'^z + {-2 + Aa)z'^, 

P-i- 

for which no recursion operator or master symmetry was known so far, and prove that this 

system admits infinitely many local higher symmetries that are constructed using a nonlocal 

two-term recursion relation rather than a recursion operator. 
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^' : Introduction 

o 

Q I Following Foursov [9J consider a Burgers-type system 

Wt = Wxx + ^wwx + (2 - 4:a)zzx, 
j^ . Zt = {I- 2a)zxx - 4:azwx + (4 - 8a)wz^ - (4 + 8a)w^z + (-2 + 4a)z^, 

where a is a real parameter. 

For a = system ([T]) is equivalent [9] to a system found by Svinolupov [33] while for a = 1 this 
system is equivalent to system (4.13) in Olver and Sokolov [26]. Finally, for a = 1/2 a recursion 
operator for ([1]) was found in [9]. 

Moreover, for any value of a the differential substitution 

Ur V 



w 



(2) 



4u 2v/^' 

the inverse of which was found in [3lj , maps the triangular system 

ut = Uxx + (1 - 2a)f^, vt = {I- 2a)vxx (3) 

into ([H). 



System ([3]) is triangular and hence C-integrable. C-integrability means here that solving (IHl) boils 
down to solving a linear heat equation for v, Vt = {I — 2a)vxx, and then solving a linear inhomoge- 
neous PDE, namely Ut = u^x + (1 ~ 2a)f ^ for a given v. More generally, C-integrability of a system 
of PDEs is tantamount to existence of a transformation that reduces this system to a linear (or at 
least triangular) one, see ^ for further details. 

In view of the above system ([T]) is C-integrable as well, and the general solution of ([T]) is obtained 
by plugging the general solution of ([3]) into ([2]). However, we are left with an open problem of whether 
the system ([1]) has infinitely many local higher symmetries for arbitrary value of a. Recall that ge- 
ometrically a higher (or generalized [25]) symmetry for a system of PDEs is essentially a solution 
of linearized version of this system; the solution in question is a vector function on the associated 
diffiety and, in general, depends on higher (i.e., of order greater than one) jets, see e.g. [3], [H], [15], 
|17j . [25] and references therein for further details. As for generalizations of this concept, see [3], 
[TT] . [13] . [H], [16], [15], [T7j and references therein for nonlocal higher symmetries, and [7], [19], [20] . 
[21j . [36] and references therein for higher conditional symmetries. 

In a somewhat more technical language, a (local) higher symmetry for ([T]) can be identified (see 
e.g. [3], [5], [11], [15], [23], [25]) with a two-component vector K = [K^^K'^Y^ whose entries K^ 
depend on x, t, w, z and on a finite number of x-derivatives of w and z up to an order k (which is in 
general different for different symmetries); here and below the superscript T indicates the transposed 
matrix. It is further required that the evolution system 

Wr = K\ Zr = K^ (4) 

is a commuting flow for ([T]), i.e., 

d'^w d'^w d'^z d'^z 



drdt dtdr drdt dtdr 
where the partial derivatives with respect to t and r are computed using ([T]) and (jlj). 

For a 7^ 0, 1/2, 1 system ([T]) was known to have [9] six higher symmetries but no recursion 
operator or master symmetry was found so far, so it was not clear whether (Q for a ^ 0, 1/2, 1 
has infinitely many (rather than just six) higher symmetries. Recall that a recursion operator is, 
in essence, an operator that maps symmetries into (new) symmetries, see e.g. [2], [5], [S], [TT], [TT], 
[2S] and references therein. In turn, a master symmetry is essentially a higher symmetry, typically a 
nonlocal one, such that the commutator of this symmetry with any given symmetry yields a (new) 
symmetry, see e.g. [2], [5], [6], [3], [21], [25] and references therein for details. 

In view of the recent results of Sanders and van der Kamp [12j who found several examples of 
two-component triangular evolution systems that possess only a finite number (greater than one) of 
local higher symmetries, it is natural to ask whether ([T]) could provide an example of a non-triangular 
system with finitely many local higher symmetries. 

In the present paper we show that this is not the case: system ([1]) has infinitely many commuting 
local higher symmetries. Quite unusual, however, is the fact that these symmetries are generated us- 
ing a nonlocal two-term recursion relation iQ rather than a recursion operator or a master symmetry, 
see Theorem [1] below for details. Moreover, we believe that it is impossible to construct a recursion 
operator of a reasonably "standard" form that would reproduce at least a part of the hierarchy from 
Theorem [1], but no proof of this claim is available so far. 



Generation of symmetries via recursion relations turns out to be of considerable interest on its 
own right. To the best of our knowledge, the first examples of this kind have appeared in [1], with the 
simplest case given by Eqs.Q, ([8]) below. However, the recursion relation ([7]) is local, and therefore 
obviously produces local symmetries; for other recursion relations in [1] locality of symmetries is also 
pretty much immediate, unlike the recursion relation ([9]) in Theorem [1] below. 

In general, if we deal with a recursion relation that involves nonlocalities, then establishing local- 
ity and commutativity of the symmetries generated using this relation is a highly nontrivial task. It 
is important to stress that the hitherto known methods for proving locality of hierarchies of symme- 
tries are based upon existence of a hereditary recursion operator or of a master symmetry, see e.g. 
[28], [31], [32], [35] and references therein, in general are not applicable in this situation. 

For the particular case of recursion relation (Q we proved locality and commutativity of the sym- 
metries Kn (IH]) using some ad hoc arguments in spirit of [11] and [18], but it would be very interesting 
to find general, more powerful methods that would not require the existence of scaling symmetry. 

It would be of interest to understand geometrical meaning of commutativity of symmetries in 
this setting. For the "standard" hierarchies this means vanishing of the Nijenhuis torsion of the 
recursion operator (see e.g. [5]) but it is not quite clear how one could generahze this to the case of 
the recursion relations. 

Finally, perhaps the most important open problem here is whether there exist any S'-integrable 
(i.e., roughly speaking, integrable via the inverse scattering transform, see 0] for details) systems 
of PDEs that have no "standard" recursion operator or master symmetry but nevertheless possess 
infinitely many symmetries generated through a recursion relation involving two or more terms. 

1 Preliminaries 

We start with recalling some basic properties of ([3]). 

First of all, for a = 1/2 this system decouples and takes the form 

Ut = Uxx, Vt = 0, (5) 

i.e., in this case we have a linear homogeneous heat equation for u, and v is simply an arbitrary 
function of x. 

On the other hand, for a ^ 1/2 upon introducing a new independent variable r = (1 — 2a)t 
instead of t in ([3]) and setting a = 1/(1 — 2a) we obtain the system that was studied, inter alia, in 
m, viz. 

Ur = aUxx + V^, Vr = V^x- (6) 

Moreover, system ([6]) (and hence ([3])) has infinitely many symmetries. Indeed, by Theorem 2.2 
of [1] the system in question has infinitely many local higher symmetries of the form 

G„ = ( ^"''" + ^« J ^ (7) 

where 6„ = bn-i — (1 — a)&n-2/2 with hi = 1, 62 = «, and 

Qn = D^Qn-l —DlQn-2 + VVn~1- (8) 
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The initial conditions for the recursion ([8]) are Qi = 0, Q2 = v"^ 
Here Ui and Vi stand for the i^^ ^-derivatives of u and v, uq e 
x-derivative (see e.g. [3], [11], [23], [25] ) 



M, f = f , and Z^a; denotes the total 



D. 



rs. 00 



9x 



i=0 



d 



d 



Ui+i- VVi+i^— 



du. 



dvi 



Recall that a function / of x, t, uq, vq, ui, vi, . 
finite number of Ui and Vi, see e.g. [3], [23] 



is said to be local if it depends only on x, t, and a 
and references therein for further details. 

Note that for a = 1 (i.e., a = 0) a recursion operator for ([6]) was found by Oevel in [21], so for 
a = we can readily find a recursion operator for ([T]) from the Oevel's recursion operator. On the 
other hand, for a = 1/2 (i.e., a = 0) Foursov [9] also found a recursion operator for ([1]). 

However, for generic a no recursion operator is known for ([6]). Moreover, for generic a the form 
of leading coefficients of higher-order symmetries of (jSD, see ([7]), appears to preclude existence of any 
"reasonable" recursion operator for ([6]), and thus for (I3|) and ([H) as well, but so far we were unable 
to prove this claim. The supposed nonexistence of recursion operator should be closely related to 
the number-theoretic aspects of the symmetry analysis for (Q, see [1], [29] and references therein. 



2 Infinitely many local higher symmetries for ([Tl) 

Using the inverse of the differential substitution ([2]) (of course, the inverse in question should be 
considered as a covering in the sense of [3], [H], [15], [16], [17] and references therein) we can obtain 
infinitely many symmetries for ([T]) with an arbitrary value of a different from 1/2 from the symme- 
tries of IQ given by ([7]), ([8]). Most importantly, all of these symmetries are local by virtue of the 
following result. 

Theorem 1 For a 7^ 1/2 system (Qj) has infinitely many commuting local higher symmetries Kn 
generated using the nonlocal two-term recursion relation 

D^ + Aw + Aw^D'^ 



Kr 



2{z,, - 2wz)D: 



Dr + 2w 



Kn-l + MKn-2 



starting from 
Here 



( 



w^ 



Kn 



Wxx + ^WWr^ 

2a -1 



+ 2zZx 



\ 



^xx \ ^W Zx 



Az{awx + (2a + l)^;^) 



2a 



2z' 



J 



M 





M 



11 



2a 



- {aDl + SawDx + 2{a{Qwx + 8w^) - (2a - l)z^ 



+ 4 {a{wxx + 8wwx) - (2a - l)zzx) D'^) 



M. 



2a 



21 



2a 



-zDx 



16a 



2a 



-wz 



2a 



2a 



{w, + Aw^)z-z^\Dz\ 



(9) 



(10) 



We defer the proof of this theorem until the next section. 

Recall (see Introduction) that locality of the symmetries Kj, j G N, means that they depend 
only on x,t,w,z and a finite number of their x-derivatives Wi and Zi {wi and Zi stand for the i^^ 
x-derivatives of w and z, and we set Wq = w and Zq = z for convenience) and do not involve any 
nonlocal quantities, cf. e.g. [3], [15], [17], [23], [25] . 

The total x-derivative D^ now takes the form 

Let / G ImDa; be a polynomial in a finite number of Wi and Zi with zero free term (let us stress 
that / is not allowed to depend explicitly on x and t). In Theorem [T] and below we make a blanket 
assumption that the result of action of D^^ on any such / again is a polynomial in a finite number 
of Wi and Zi with zero free term, i.e., we always set the integration constant to zero. Note that an 
alternative (and more general) approach is to define the operator D^^ in the fashion described in 

[ID], [22], [SO]. 

With the above assumption in mind we readily find that the first nontrivial symmetry generated 
via ([9]), namely, K^, has the form K^ = {K^, K^)'^, where 



2 



, (« + l) 12(a + l) ^ 12(a + l) 

iio = Wxxx WWxx + OZZxx Ul 

3 2a - 1 2a - 1 2a - 1 ^ 

+ 3zl — Qz^Wx w'^Wx + I2wzzx, 

2a — 1 

r^2 6a ^ ^ 12(4a + l) ^,^ ^ 2^ 

K^ = Zxxx + 7: tZWxx + OWZxx + QWxZx H r ; WZWx + 6{2w - Z )Zx 

Za — i 2a — i 

, 24(2a + l) 3 ^^ 3 
2a -1 

We see that the symmetries Kj, j = 1,2,3 are independent of x and t. Moreover, it is easily seen 
that by virtue of the above definition of action of D~^ the symmetries Kj in fact do not depend 
explicitly on x and t for all j G N. 

3 Proof of Theorem [T] 

The recursion relation ([8]) can be rewritten in the terms of symmetries Gi from ([7]) as 

G„ = D.(G„_0 - ( ^^ ~^^'^ ° \ Dl{G^_,) +(00) Gn-.. (11) 

and the recursion relation ([9]) readily follows from ( TTT1) and ([2]). 

Now we have to show that the symmetries Kj generated via the recursion relation ([9]) with the 
initial data (fTOj) are local in the sense of the preceding section for all j = 3,4, . . .. This obviously 
boils down to proving that they do not involve nonlocalities like y = D^^{w). 

To this end we shall use induction on n starting from n = 1. In view of the specific form of the re- 
cursion relation (^ we only have to show that once Kj_i = {Kj_^, Kj_^)'^ and Kj_2 = (-^j-2' ^j-2)^ 



are such that Kj_i, Kj_2 G Im-Dj; (and hence Kj is local) then we have Kj G linDx- This is obviously 
true for j = 3, so we only need to establish validity of the induction step. This is done using the 
following result. 

Lemma 1 Fix j > 3 and assume that Ki, i = j — 2,j — 1, are local, and Kj_i, Kj_2 G IraDx- Then 
Kj defined via ^ is local too, and we have Kj G ImD^. 

Proof of the lemma. As Kj_i, Kj_2 ^ Ihi-D^: by assumption, locality of Kj is immediately inferred 
from (ini), so it remains to prove that Kj G IraDx- 

Now, the condition Kj G ImDj. can be restated as follows: pj G ImD^, where pj = p'q[Kj], Pq = w 
is a conserved density for ([1]), and f'[K] stands (cf. e.g. |2], [S], [21]) for the directional derivative 
(also known as linearization, see e.g. [3], [I5], [H]) of / along K = {K^,K'^)^: 



i=0 ^ ^ 



Recall (see e.g. [3], [TT], j23], [25]) that a local function p is said to be a local conserved density ioT. ([T]) if 
Dt{p) G ImDx. A local conserved density p is said to be nontrivial if p ^ lavDx, and trivial otherwise. 

Next, as Kj is a symmetry of ([1]) by construction, the quantity pj is a local conserved density for 
(HD, see e.g. [2], [11], [3], [25]. 

However, for a 7^ 1, 1/2 the quantity po is the only nontrivial local conserved density for ([1]) mod- 
ulo the terms from \yivDx- Indeed, by virtue of the results of [23] (see also Theorem 5-1 of [8], and cf. 
[27] and references therein) any nontrivial local conserved density for ([T]) depends (again modulo the 
terms from ImD^) only on x, t, w, z, w^, z^, w^x, Zxx, and the direct search for all conserved densities 
of this form proves our claim. 

Hence the most general local conserved density for ([T]) has the form cpo + p, where p G ImDx and 
c is a constant. In particular, if Kj is local then we have 

Pj = CjPo + Pj, (12) 

where pj G ImDx and Cj are constants. 

In analogy with Krasil'shchik [18], let us show that if Kj is local then Cj = 0, and thus pj G ImDx. 
First of all, note that system ([1]) has [9j a scaling symmetrjo 

(w 

and it is readily seen that all Kj constructed using ([9]) are 5-homogeneous: we have Ls{Kj) = 
[S, Kj] = jKj, j = 1, 2, . . .. Here Lq denotes the Lie derivative along Q, see e.g. [2], [25], [28], [32] 
for details, and [■, ■] is the standard commutator of symmetries (see e.g. [2], [3], [5], [TT], [17], [23], [25]): 

[F,G] = G'[F]-F'[G]. (13) 



^Note that S, Ki, and 1^2 exhaust all linearly independent Lie point (or, more precisely, higher symmetries that 
are equivalent to the Lie point ones in the sense of [TT1[25]) symmetries of ([T]). 



Hence Ls{pj) = p'j[S] = (j + l)pj + rjj, where rjj G ImDx. On the other hand, 

Lsicjpo + Pj) = Cjpo + (j, 

where Q G IraD^. Therefore, if we act by Ls on the left- and right-hand side of ( TT2l) and equate the 
resulting expressions, we obtain 

jcjpo + Oj = 0, (14) 

where 6j G ImD^. 

As po ^ Iih-Dx and j 7^ by assumption, Eq.( fT4l) can hold only if Cj = 0. Hence pj = pj G ImD^., 
and the lemma is proved. D 

Now that we have proved locality of Kj for all j G N, we can easily prove commutativity of Kj 
with respect to the bracket (TT3l) : [Ki, Kj] = 0, i,j = 1, 2, 3, . . .. 

First of all, in analogy with the reasoning presented in Chapter 4 of [11] for scalar evolution 
equations (cf. also [3], [I8j, [I5]) it can be shown that any S-homogeneous x, t-independent local 
higher symmetry G of ([I]) of order j > 1 is of the form 

ajWj + G] 

G=\ .'^ \, (15) 

(3jZj + Gj 

where aj, Pj are constants, and Gj, G'j are polynomials in w, z, wi, zi, . . . , Wj-i, Zj_i, and these poly- 
nomials have no free terms and no linear terms. 

Let Ck be the space of S-homogeneous x, t-independent local higher symmetries of ([I]) of order 
no greater than k. By the above, there exists a basis in Ck that consists of symmetries of the form 
flTSj) for j = 1,2, ... ,k (note that in principle this basis may contain more than one symmetry of 
given order j). 

Clearly, Ki G Ci and [Ki,Kj] G Ci+j. But using ( TT3l) and ( TT5l) we readily see that the commu- 
tator [Ki, Kj] contains no linear terms, and hence (cf. e.g. [H]) this commutator can belong to Ci+j 
only if [Ki, Kj] = 0. Thus, [Ki, Kj] = for all i,j G N, and this completes the proof of Theorem [H 
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